Inter mittency analysis in diffractive DIS processes 

Zhang Yang 

Institut fiir theoretische Physik, FU Berlin, Arnimallee 14, 14195 Berlin, Germany 

It is pointed out in this talk that "the colorless objects" in diffractive scattering in the small-ss region can 
be probed by measuring the scaled factorial moments of fmal-state-hadrons and the dependence of their scaling 
behavior upon the diffractive kinematic variables. The Monte Carlo implementation of RAPGAP and JETSET 
are discussed as illustrative examples. The results of these model-calculations show in particular that inclusion of 
the contributions from the photon gluon fusion processes can considerably enlarge the power of the scaled factorial 
moments. The significant IB-dependence of intermittency index in small- xb region in this MC analyses shows 
the hadronic final states depend very much on the lifetime of the exchanged colorless object Cq. The dependence 
of the intermittency index upon ip, Q 2 and /3 are also discussed. 



Recent years have witnessed a remarkably in- 
tense experimental and theoretical activity in 
research of the dynamics of diffractive deep- 
inelastic electron-proton (ep) scattering (DIS) in 
low Bjorken-x kinematic range. The measure- 
ments at HERA have demonstrated the existence 
of a distinct class of events in which there is no 
hadronic energy flow in an interval of pseudo- 
rapidity adjacent to the proton beam directionjp. 
Our present understanding of DIS could be inade- 
quate at Iow-xb because additions to the leading 
order QCD-based partonic picture are likely to 
be substantial. A natural interpretation of these 
so-called "large rapidity gap" events is based on 
the hypothesis that the deep-inelastic scattering 
process involves the interaction of the virtual bo- 
son probe with a colorless component of the pro- 
ton. Hence there is no chromodynamic radia- 
tion in final state immediately adjacent to the 
direction of the scattered proton or any proton 
remnant. What is this colorless component orig- 
inating from the nucleon? The large rapidity 
gap events discovered in deep inelastic scatter- 
ing at HERAjl[ are usually interpreted in terms 
of pomeron-exchanged model Although this 
seems to work reasonably well phenomenologi- 
cally, there is yet no satisfactory understanding of 
the pomeron structure and its interactions mech- 
anism. In this respect, it is helpful to probe the 
properties of such objects in the non-traditional 
aspects, in addition to the traditional measures 
such as rapidity gap distribution, structure func- 



tion and the averaged cross section and so onjlj. 
In the present talk we wish to point out that use- 
ful information about the exchanged colorless ob- 
ject in diffractive lepton- nucleon scattering pro- 
cess can be extracted by studying scaling behav- 
ior and fractality (intermittency) of the final state 
of the colorless component in proton excited by 
the virtual boson probe. 

The manifestation of fractality (intermittency) 
in high energy multiparticle production process is 
the anomalous scaling |^ 

f Ax\ 4 ' q 

F q (Sx) = F q (Ax) ( — J , as M -> oo, Sx -> 0(1) 

of g-order factorial moments (FM's) F q , defined 
as 

Fg (6x) = ±Y Mn rn -l) (n m - q + l)) 

where x is some phase space variable, e.g. 
(pseudo-)rapidity, the scale Sx = Ax/M is the 
bin width for a Af-partition of the region Ax 
in consideration, n rn is the multiplicity in the 
mth bin. Since the factorial moments can re- 
move the statistic noise around the probability 
and associated directly with the scaled moments 
of probabilities ||, the scaling exponent <fr q in 
Eq.( [l]), called intermittency index, characteriz- 
ing the strength of dynamical fluctuation, is con- 
nected ji) with the anomalous fractal dimension 
d q of rank q of the spatial-temporal evolution of 
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high energy collisions, d q — 4> q /(q — 1). A possi- 
ble cause leading to the power-law of FM's of final 
particles in high-energy collisions is that the emit- 
ting source of final hadrons is self-similar frac- 
tal||. Furthermore, DIS experiments and the em- 
pirical analyses 0,0] show that the gluon-density 
in the nucleon in the low-irs kinematical region 
is much higher than those for quarks/ antiquarks, 
and it is increasing with decreasing i£_b(7]]. In 
this soft gluon system where the gluons inter- 
act with each other in complicated dynamic pro- 
cesses, it has been proposed^ that the gluons 
may 'self-organize' into the clusters in the dissi- 
pative gluon system by self-organized criticality 
(SOC)||, and the colorless component in pro- 
ton can be regarded as color singlet gluon clus- 
ters The spatial- temporal structure of the 
SOC-cluster (or BTW-cluster) is self-similar frac- 
tal||. So it is feasible to study the fractal struc- 
ture of the colorless component of the proton by 
measuring the anomalous scaling behaviors of fac- 
torial moments of the final state particles origi- 
nating from the scattering of virtual photon and 
colorless object in the diffractive lepton- nucleon 
scattering. 

The xs-dependence of scaling behaviors of 
FM's, and in particular that of intermittency in- 
dex 4> q (Q 2 , xb', xjp) for fixed Q 2 (and a;p), plays 
a distinguished role in such studies. This is be- 
cause, when the virtual photon 7* originating 
from the incident lepton is absorbed by the nu- 
cleon, its energy- momentum q = (q° ,q) is in fact 
absorbed by a virtual colorless component of the 
nucleon in the lepton-nucleon diffractive scatter- 
ing. In a fast moving reference frame, for example 
the lepton-nucleon center-of-mass frame, where 
the nucleon's momentum P is large in high-energy 
collisions, the time interval Ti n t in which the ab- 
sorption process takes place (it is known as the 
lepton-nucleon interaction/collision time) can be 
estimated by making use of the uncertainty prin- 
ciple. In fact, by calculating l/q° in this reference 
frame we obtain 1 11 1: 



This means, for given P and Q 2 = —q 2 , the in- 
teraction time ri n t is directly proportional to is 



in the Iow-xb {xb "C 1) kinematic range. In 
other words, is a measure of the time- interval 
in which the absorption of 7* by the space-like 
virtual colorless object takes place. Hence, by 
studying the ^-dependence of the intermittency 
index <fi q (Q 2 , xb', sp), we are not merely probing 
the anomalous scaling behaviors of the collision 
process between 7* and the colorless object which 
we hereafter call Cq . Since this hadronization pro- 
cess of the virtual colorless object Cq is initiated 
by the interaction with 7*, we are also examing 
whether/how the hadronization process changes 
with the interaction time Tint. This question is 
of considerable interest, because a virtual photon 
7* can (logically) only be absorbed by virtual sys- 
tems (cg's) whose lifetimes (r*'s) are longer than 
interactive time Tint (i.e. Tint < r c ). That is, the 
average lifetime (r c ) of the Cq's, which can ab- 
sorb a 7* associated with interaction-time Tint , is 
a function of Tj nt . Hence, from the xs-dependence 
of the scaling behavior of FM's, in particular from 
that of the corresponding 4> q (xB, Q 2 ] xp)'s, we 
can in principle find out whether/how the dy- 
namic fluctuation and the fractal structure in the 
spatial-temporal evolution originating from the 
exchanged colorless object Cq depends on its av- 
erage lifetime (r c ) of the Cq's. 

The Q 2 -dependence of the scaling behavior of 
the FM's is also of considerable interest. This 
is because, in photon-proton scattering experi- 
ments, not only those with real (Q 2 = 0) photons 
but also those with space-like (Q 2 > 0) photons 
where Q 2 is not too large (< lGeV 2 /c 2 , say) have 
very much in common with hadron-hadron colli- 
sions. Having in mind that the index of inter- 
mittency for hadron-hadron scattering is smaller 
than that for electron-positron annihilation pro- 
cesses Q, we are led to the following questions: Do 
we expect to see a stronger Q 2 -dependence when 
we increase Q 2 from zero to 10 or 100 GeV 2 /c 2 , 
say? Is this also a way to see whether space-like 
photons at large Q 2 " behave like hadrons" in such 
interactions? 

While waiting for data to perform the above- 
mentioned analysis, let us now consider the fol- 
lowing two phenomenological models as illustra- 
tive examples: 

(A). If the colorless object (cq) is a quark- 
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antiquark pair (formed by interacting gluons) 
which exists in the time-interval when the virtual 
photon 7* is absorbed by the object, we shall see 
the following: Especially when Q 2 is sufficiently 
large, the incoming 7* (the transverse dimension 
is expected to be proportional to l/Q 2 ) will hit 
the quark (g) or the antiquark (q) and make them 
fly apart symmetrically with respect to the cen- 
ter of mass the j*qq system — similar to the qq- 
pair produced in e + e~-collisions (with respect to 
the center of mass the qq system). That is, in 
this case, the final-state-hadrons of an event are 
fragmentation- products of the quark and/or the 
antiquark, and hence they are expected to show 
characteristic features similar as those observed in 
the reaction e + e~ — > hadrons at the same c.m.s. 
energy. It is interesting to see that the very recent 
inclusive measurements performed at HERAjl]] 
in diffractive electron-proton scattering show the 
following: Both the scaled longitudinal momen- 
tum (xp) distribution and the transverse momen- 
tum (p±) distribution are strikingly symmetric 
with respect to the center-of-mass of the photon 
and the struck colorless object; and the general 
features of these distributions are very much the 
same as those observed in electron-position colli- 
sion processes. These facts strongly suggest that 
a more detailed comparison between these two 
collision processes would be useful. 

For this purpose, we made use of the Monte 
Carlo (MC) program JETSET||. We gener- 
ated 50,000 MC events, and calculated the second 
normalized factorial moment in 3-dimensional 
(i~),p±,<fi) phase space at the given qq-c.m.s en- 
ergy y/s, where the pseudorapidity i], transverse 
momentum p± and the azimuthal angle <j) are 
defined with respect to the sphericity axis of 
the event. The usual cumulative variables X 
translated from x — (rj,px,(f)), i.e. X(x) = 
j x p(x)dx/ J x max p(x)dx, were used to rule out 
the enhancement of F q from a non-uniform inclu- 



sive spectrum p(x) of the final hadrons|12|. The 
obtained results have been collected in Fig. la in 
double logarithmic F2 vs M plots for different 
y/s (or Mx , which is the invariant mass of the 
7*Cg system in a corresponding diffractive lepton- 
nucleon scattering event). It is clear that the 
higher the invariant mass of 7*Cq is, the larger 



the power of FM's, i.e. the stronger the dy- 
namic fluctuation is. In the very low invariant 
mass (y/s = 4.5GeV, say), the powers of FM's 
become less than 0, which can be referred to the 
constraint of the momentum conservation in the 
hadronization process p4). 





Vs = 4.5 GeV 

" o 2 =-0.234±0.002 


Vs = 7.5 GeV 

" ffl 2 =-0.110±0.002 


= 10.5 GeV 

" o 2 =-0.065±0.002 


vs = 13.5 GeV 
" o 2 =-0.032±0.002 




■ 








Vs = 16.5 GeV 

. ...^ 

. 1 |>. ! =0.006±0.002 


Vs = 19.5 GeV 
, •— 

_ e 2 =0.025±0.002 


Vs = 22.5 GeV 
. o 2 =0.067±0.002 


vs = 25.5 GeV 
. o 2 =0.096±0.003 


vs = 28.5 GeV 
. o 2 =0.125±D.003 


2 


2 


2 

(a) 


2 


2 

(Ln M)/3 


.0 < M x < 3 (GeV) 
N„ venl =85647 


.3 < M x < 6(GeV) 
N even| =1 76652 


.6 < M x < 9 (GeV) 
N oven| =1 05763 


3 <M X <12(GeV) 
N event =5461 8 


.12 <M X < 15 (GeV) 
N evl , nt =29074 


o 2 =0.174± 0.002 


o 2 =0.263± 0.002 


o 2 =0.34B± 0.002 


ij> 2 =0.399± 0.002 


4i 2 =0.431± 0.002 

>* 


77 * ' 


: 

y • 




7 




.15 <M„ < 18 (GeV) 


.18 < M K <21 (GeV) 


.21 <U X < 24 (GeV) 


.24 < M x < 27 (GeV) 


27 < M x < 30 (GeV) 


N ovenl =165B8 

: ./ 


N evenl =10094 

V 


N evenl =6656 


N elrent =41B7 , 

/ 

7 


N elr( . nt =2840 * I 

/ 

/ 


7 

' / >2=0-468+ 0-002 


' f* <(i 2 =0.486± 0.002 


'/* i|> 2 =0.485± 0.003 


7 * *2=0-535± 0.003 


■ /• $ 2 =o. 540 0.003 


2 


2 


2 

(b) 


2 


2 
(Ln M)/3 



Figure 1. The scaled factorial moments 

F2 versus the number M of subintervals of 3- 
dimensional (f],p±, <j)) phase space in log-log plot, 
and the second-order intermittency index 02 in 
corresponding sample set. (a). The MC result 
of JETSET 7.4||!] in different cms energy y/s, 
50000 events are generated in each sample set; 
(b). That of RAPGAP||| in corresponding in- 
variable mass Mx interval with N evQn t events in 
each subsample. 

(B). Based on pomeron-exchanged modcl^), 
several Monte Carlo generators, such as 
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POMPYT@, RAPGAPjT|l have been set up. 
These models have been used to describe quite 
well the HERA data in the global and aver- 
aged features, such as the rapidity gap distribu- 
tion, the diffraction and total cross section and 
so onjjj. The theoretic models are usually con- 
fronted with incorrigible discrepancy, when the 
data about the locally nonstatistic fluctuation in 
small phase space are involved in the compari- 
son with themQ. No evidence has shown that 
this kind of locally dynamic fluctuation could 
be certainly referred only to the hadronization 
process but have nothing to do with the initial 
stage of high energy collisions. In this respect, 
it is also relevant to see whether the features 
of scale invariance and fractality in small phase 
space of the lepton-nucleon diffractive processes, 
specially their dependence upon the diffractive 
variables can be compatible with the pomeron 
type of model. In the following we take RAP- 
GAP as an example, in which the virtual photon 
(7*) interacts directly with a parton constituent 
of the pomeron either in lowest order (Fig. 2a) or 
in 0(a oul a s ) order — via the photon gluon fu- 
sion (Fig. 2b). In both cases a color octet rem- 
nant is left at low pt with respect to the pomeron 
and hence also to the proton. The higher order 
gluon emission is simulated with the Color Dipole 
Model (ARIADNE) JTtJ and the hadronisation is 
performed using the JETSET jl3]]. The pomeron 
flux factor fpp(t, xp) and pomeron structure 
function G(/3) are taken respectively as 

/pp(^p) = %%- 2Qp(t) , (4) 
given by Berger et al. and Streng|^], and 
pG{fi) = 60(1 - (5) 

when the pomeron is made of 2 'unrealisti- 
cally hard' gluons as suggested by Ingelman and 
SchleinQ. Here the kinematic variables t, xp 
and j3 are defined as t = (P — P') 2 , xp = 
q-(P- P')/(q ■ P), and (3 = -q 2 /[2q ■ (P - P% 
where P' is the 4- momentum of the final state col- 
orless remnant of nucleoli. The pomeron Regge 
trajectory is given by ap(t) = 1+e+a't, e w 0.085 
and the slope a' = 0.25 obtained from a fit to 
datag. 




pp. pp. 



(a) (b) 

Figure 2. The basic processes included in the 
RAPGAP [fl6[ implementation for inelastic lepton 
scattering on a pomeron: (a) the lowest order pro- 
cess for hard parton level; (b) the 0(a cm a s ) order 
process for photon gluon fusion. 



We generated 500,000 RAPGAP events, and 
divided the whole sample into 10 subsamples ac- 
cording to the invariable mass Mj of 7*Cq in the 
MC event. The scaling behaviors of FM's for dif- 
ferent Mx intervals are shown in Fig. lb. The de- 
pendence of scaling behaviors of FM's upon Mx 
is similar with the result of JETSET in Fig. la, 
i.e. the powers increase for increasing scattering 
energy of 7*cJ. But it is noticeable that the inter- 
mittency index fa for a given Mx interval is much 
larger in RAPGAP than that in JETSET. Hav- 
ing in mind that in Fig. la the colorless object 
is considered as a quark-antiquark pair formed 
by gluons and the Feymann graph of the diffrac- 
tive process in this aspect is just same as shown 
in Fig. 2a, the difference between cases (A) and 
(B) is that the higher-order photon-parton inter- 
action has been taken into account in RAPGAP 
(see Fig. 2b). It is understandable since the 
branch number of the parton cascading process 
in parton shower level is larger when the pho- 
ton gluon fusion is involved in Fig. 2b, while it 
is believed generally thatH power-law behaviors 
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Fig. 3 

Figure 3. The dependence of 2-order inter- 
mittency index <p 2 in RAPGAP jl6| Monte Carlo 
implementation upon diffractive kinematic vari- 
ables: (a) xb, (b) Q 2 , (c) xp and (d) (3. 



in the color-string type of models are referred in 
large part to the randomly cascading process of 
parton energy in perturbative phase, so the frac- 
tal in the evolution processes with longer cascade 
branch would be stronger. 

In order to see the dependence of dynamic fluc- 
tuation upon the other diffractive variables as we 
argued above, we divided the RAPGAP MC sam- 
ple into 10 subsamples according to and Q 2 . 
The results of scaling behaviors of FM's are shown 
in Figs. 3a and b respectively. This example ex- 
plicitly shows how the proposed method can be 
used to analyze the feature of the evolution pre- 
cess of the colorless object in diffractive lepton- 
nucleon scattering processes. The significant in- 
dependence of the intermittency index observed 
in this example not only shows that the hadronic 
final states depend very much on the average life- 



time of the exchanged colorless object Cq if Cq 
can be indeed considered as a pomeron as sim- 
ulated in RAPGAP, but also the following: The 
longer the average lifetime of Cq is, the stronger 
the fractal structure of the space-time evolution 
processes originating from the 7*Cq's interaction 
will be. The smaller the Q 2 becomes, i.e. the 
larger transverse dimension of virtual 7* is, the 
more 7* "behaves like a hadron" . 

Since ip may be interpreted as the fraction of 
the 4-momentum of the proton carried by the ex- 
changed colorless objects, and (3 as the fraction 
of the 4-momenta of the exchanged objects car- 
ried by the parton interacting with the virtual 
boson, it meaningful to calculate the dependence 
of scaling behavior upon xp and /3 (Fig. 3c and 
d respectively). It seems from RAPGAP that 
only when the momentum of colorless object is 
large enough, can scaling behavior of the final 
state originating from 7*Cg be significant; and 
it is also clear from Fig. 3d that the dynamic 
fluctuation in the diffractive scattering don't in- 
crease monotonously for increasing the fraction of 
energy of stricken parton in the pomeron. 

In conclusion, we have shown in this note 
that the colorless object Cq in diffractive lepton- 
nucleon scattering can be probed in a model- 
independent way, i.e. by performing the scaled 
factorial moment analyses for the hadronic final 
states originating from the scattering of virtual 
boson 7* and the Cq on an event-by-event bases. 
The Monte Carlo implementation of RAPGAP 
and JETSET have been discussed as illustrative 
examples. The results of these model-calculations 
show in particular that inclusion of the contribu- 
tions from the photon gluon fusion processes can 
considerably enlarge the power of the scaled fac- 
torial moments. The significant a^-dependence 
of intermittency index in small- xq region in this 
MC analyses shows the hadronic final states de- 
pend very much on the lifetime of the exchanged 
colorless object cj$, i.e. the longer the lifetime 
of Cp is, the stronger the fractal structure of the 
space-time evolution processes originating from 
the 7*Cq's interaction will be. 
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